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G. J. Piermarini* and S. Block 
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An X-ray diffraction study of a barium borosilicate glass consisting of 24 mole percent 
barium oxide, 40 mole percent boric oxide, and 36 mole percent silicon dioxide has been per- 
formed. Resulting atomic radial distribution functions give the following average inter- 
atomic distances: Si-O, 1.6 A; Ba-O, 2.8 A; Ba-Ba, 4.7 A; and Ba-Ba, 6.8 A. From the 4.7 
A Ba-Ba separation a Ba-O-Ba bond angle of about 115° has been calculated. The observed 
average barium separations are in partial agreement with that predicted by Levin and Block 
on the basis of a structural interpretation of immiscibility data. A proposed coordination 
change by Levin and Block for the barium atoms in the system has been confirmed, but the 
details of the coordination change mechanism have not. 

Combining the results of the radial distribution study and immiscibility data on the 
barium borosilicate modifier-rich Liquid at maximum barium oxide content has indicated 
that approximately 16.75 mole percenl barium oxide is involved in the 4.7 A separation and 
8.25 mole percent is associated with the 6.8 A separation. 

A mechanism which allows the composition of the modifier-rich liquids in I he ternary 
system to be calculated has been presented. The calculated composition has been found to 
agree well with the experimental value. 



1. Introduction 

A comprehensive structural interpretation of the 
immiscibility of alkali and alkaline-earth metal 
oxides in borate, silicate, and borosilicate systems 
has been proposed by Levin and Block [l]. 1 Two 
models representing short-range structures, orig- 
inated by Warren and Pincus [2], were adopted to 
enable the composition of the immiscibility limit 
to be calculated. The adopted models have been 
described previously [1,2] and are shown in figures 
1 and 2. 

For the binary barium borate system at the limit 
of maximum miscibility of barium oxide in the modi- 
fier-rich liquid, Levin and Block proposed that the 
average separation of barium atoms is 6.67 A 
resulting from type B coordination. Subsequent 
radial distribution analysis of the appropriate 
barium borate glass by Bienenstock, Posner, and 
Block [3] demonstrated that the proposed model 
was essentially correct. The average separation of 
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Figure 1. Schematic representation of type A coordination. 

The network former is in either tetrahedral or triangular configuration. The 
calculated barium-barium nearest neighbor separation (dash line) is 5.50 A. 




Figure 2. Schematic representation of type B coordination. 

The network former is in tetrahedral configuration. The network former and 
barium polyhedra share edges. The calculated barium-barium nearest neighbor 
separation (dash line) is 6.67 A. 



barium atoms was found to be 6.8 A. This conclu- 
sion was based on the assumption that the structure 
of the modifier-rich liquid did not differ appreciably 
from the structure of the corresponding glassy state 
obtained by rapid quenching [4]. 

For the ternary system, barium borosilicate 
(part A, fig. 3), it was postulated that the configura- 
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Figure 3. Immiscibility regions {weight percent basis) in 
systems of type RO— B 2 03-Si02 showing form of modifier-rich 
boundary curves. 

Numbers in parentheses, adjacent to intersection of solubility curves and 
tangents, indicate maximum solubility of modifier. To explain the form of these 
boundary curves, a change in coordination type for the modifier cation has been 
proposed by Levin and Block. Reproduced from E. M. Levin and S. Block, 
J. Am. Ceram. Soc. 40, 95 (1957). 



tion of the barium atoms in the modifier-rich liquids, 
at the limit of maximum miscibility of modifier 
oxide, changed progressively upon addition of 
silicon dioxide component from type B coordination, 
characteristic of the barium borate system, to type 
A coordination characteristic of the barium boro- 
silicate system at maximum barium oxide content 
(see parenthetical number in part (A) of fig. 3 and 
related legend). Accompanied with the coordina- 
tion change, a decrease in the average separation 
of barium atoms from 6.67 A (type B coordination) 
to 5.50 A (type A coordination) was postulated. 
The primary purpose for this research was to 
compare the average barium interatomic separation 
in barium borate and barium borosilicate systems 
at maximum miscibility of barium oxide in the 
respective modifier-rich liquids in order to deter- 
mine the mechanism accounting for the observed 
increase in the immiscible range of barium boro- 
silicate. This paper describes the determination 
of average interatomic separations by means of a 
Fourier analysis of X-ray diffraction patterns of the 
appropriate barium borosilicate glass. It also pre- 
sents a qualitative interpretation of the observed 
interatomic distances. With respect to the average 
barium separation, if the coordination change 
(B to A) postulated by Levin and Block is a physical 
reality, then the resulting atomic radial distribution 
curve should show the presence of a large maximum 
in the 5.5 A region and the absence of a peak of 
comparable magnitude in the 6.8 A range. 



2. Experimental Procedures 

The experimental procedures were essentially 
the same as those employed previously in the vitreous 
barium borate study [3] . Measurements were made 
with a Norelco X-ray diffractometer. Intensities 
were determined with a Norelco counting rate com- 
puter adjusted so that the probable intensity error 
did not exceed 2.4 percent. Monochromatization 
of the scattered radiation was achieved by means of 
the balanced filter technique [5]. Intensity cor- 
rection due to air scatter was avoided by evacuating 
the volume defined by the scatter shield surrounding 
the specimen. The intensity data were auto- 
matically determined and recorded on a Brown 
strip-chart recorder at 0.05° intervals in 20. Using 
Cu radiation, measurements were taken in the range 
2° to 70° 20. With Mo radiation, the intensity 
measurements were made in the range 20° to 60° 
20. Periodic checks were made to guard against 
significant variations in incident intensity. Further 
detail in the experimental procedure is given in 
reference 3. 

The specimen, containing 24 mole percent barium 
oxide, 36 mole percent silicon dioxide, and 40 mole 
percent boric oxide, was obtained from G. W. Cleek 
and E. H. Hamilton of the National Bureau of 
Standards staff. The measured density of this 
specimen, determined by the method described by 
Glaze, Young, and Finn [6], is 3.35 g cm -3 . The 
normalized experimental intensities (A) corrected 
for polarization are shown in figure 4. 

The coherent scattering factors, /, for barium, 
boron, silicon, and oxygen, in addition to the in- 
coherent scattering factors for silicon, boron, and 
oxygen, were obtained from standard references [7]. 
The incoherent scattering factor for barium was 
calculated using [8] 

l lne =Z-f/Z. 
The incoherent scatter was corrected with the 
Breit-Dirac factor [8] . The total independent scatter 
curve (B) in figure 4 is the sum of the coherent (C) 
and incoherent (D) scatter. Unlike the technique 
employed in the barium borate study [3], the experi- 
mental intensities corrected for polarization were 
put on an absolute scale (electron units per unit 
of composition) by adjusting the ordinates of the 
experimental intensity curve so that the experimental 
curve (A) matched the total independent scattering 
curve (B) at large s values. This procedure was 
followed because it appeared less likely to produce 
high frequency ripples in the distribution curve 
(assuming that systematic errors in intensity meas- 
urements and scattering curves have not been 
completely eliminated) which would interfere with 
the determination of peak positions. The method 
used here tends to push the error residue toward 
low s values thereby producing longer wavelength 
ripples which interfere less with the determination 
of peak positions. That high frequency ripples 
lead to inaccuracies in the distribution curve seems 
to be indicated in the barium borate study [3]. 
The normalized intensity curve (A) corrected for 
polarization is shown in figure 4. 
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Figure 4. Normalized experimental X-ray intensities (A) and 
independent scattering curves (B,C,D) for vitreous BaO- 
B 2 3 -Si0 2 . 

Following essentially the same calculations em- 
ployed in the previous barium borate study [3], 
the atomic radial distribution functions were de- 
termined. The weighted distribution is usually 
given by the relation: 



S#»4iir*p w (r) = IX™47r/' 2 p + 



2r r 

TTJS r 



si(s) sin sr ds, 



where 2 indicates summation over the molecular 
composition. K m is the effective number of electrons 
in atom m, and is defined as K m =f m /f e . 4irr 2 p m (;r)dr 
is the number of atoms, each multiplied by its 
effective number of electrons, between distances r 
and r-{-dr from atom m. p is the average number 
of electrons per unit volume, s equals (4x sin 6) /A, 
where 6 is the Bragg angle of diffraction and X 
is the wavelength of the radiation. i(s) is the experi- 
mental amplitude function and is equal to (I eu — 
*Zjm)lj 2 e, where I eu is the corrected experimental 
intensity of unmodified scattering in electron units 
per unit composition, f m is the atomic scattering 
factor for atom m, and f e is the average/ per electron. 
That is, f e equals S/ m /SZ m , where Z m is the atomic 
number of atom m. The reader is referred to refer- 
ence 9 for further mathematical detail related to 
the above general equation. 

The determination of the radial distribution 
function comprises two main steps: First, the 
numerical evaluation of the si(s) function from 
experimental scattering data, and second, the 
evaluation of the summation 

2r Smax 

— 2 s i( s ) sm srAs. 

7T Sn.in 



These summations were performed on a high speed 
digital computer for 0.1 A intervals in r from r=0 
to r=9.9 A over specified ranges in s, giving directly 
the differential radial distribution function, 

J2K m 4irr 2 [p fri ^)-p ], 

from which, by adding the average radial atomic 
density, 2K m 4irr 2 p , the atomic radial distribution as 
a function of r is obtained. 

3. Results 

A fit at high s values of the corrected experimental 
intensity (A), derived from both the molybdenum 
and copper radiations, to the total independent 
scatter curve (B) is shown in figure 4. An atomic 
radial distribution function derived from the experi- 
mental intensity data is shown in figure 5. Several 
techniques were employed in determining the 
radial distribution function. For example, the 
series was terminated at various s values, the 
experimental scatter curve (A) was extrapolated 
to zero at s=0, the CuK a and MoK« intensity 
curves were fitted together in several ways, a 
number of normalization trials were made, and 
various temperature factors were applied to the si(s) 
functions. Application of these techniques produced 
no significant change in the peak positions in the 
distribution curve shown in figure 5. However, the 
peak areas were affected. A temperature factor of 
e -o.o2s 2 was a ppli e d to the experimental amplitude 
function used in determining the distribution 
curve shown in figure 5. The intensity data covered 
the range 0.4 to 8.8 A -1 in s. 
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Figure 5. Atomic radial distribution function of Ba02-B20i 
Si02 glass. 
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J? or the purpose of this study the more important 
features of the distribution function shown in figure 5 
are the presence of two unexpected major peaks, one 
at 4.7 A, the other at 6.8 A, and the absence of an 
expected major peak in the 5.5 A region. The two 
observed peaks are attributed primarily to barium- 
barium interatomic separations, because the barium 
atom is the only one present in the glass having a 
large enough scattering factor to produce peaks of 
this magnitude. 

The 4.7 A barium -barium separation is associated 
with a modified type A coordination wherein the 
barium-oxygen-barium bond angle is approximately 
115, not 180° as postulated by Levin and Block. A 
schematic representation of modified type A coordi- 
nation is shown in figure 6. If it is assumed that 
the barium atoms form, in general, a cubic array 
with an average edge of 4.7 A, then the face diagonal 
of this cube has a length of 6.65 A. Although the 
radial distribution function shows high atomic den- 
sity in the region of 6.65 A, the maximum is at 6.8 A. 
This peak is primarily due to barium atoms in type 
B coordination, and the second nearest neighbor 
barium -barium separation arising from modified type 
A coordination. A similar barium separation (6.8 A) 
was found for type B coordination in the barium 
borate study by Bienenstock, Posner, and Block [3]. 

Further analysis of the distribution function indi- 
cates that the 1.6 A peak represents primarily the 
average silicon-oxygen nearest neighbor distance. 
The calculated distance is 1.62 A. By a similar 
radial-distribution study, Warren et al, obtained 
1.62 A for the silicon-oxygen nearest neighbor dis- 
tance in silica glass [9]. The 2.8 A peak represents 
the barium-oxygen nearest neighbor distance. The 
calculated separation is 2.8 A. 

Coordination numbers calculated from the ob- 
served areas of maxima are rather doubtful due to 
systematic errors and the overlap of interatomic sepa- 
rations inherent in the technique. In addition, for 
complex glasses such as barium borosilicate, the cal- 
culation of coordination numbers is complicated. 

Consider the 1.6 A peak. There are three inter- 
atomic separations contributing to the area of this 
peak. One is the previously mentioned 1.6 A silicon- 




Figure 6. Schematic representation of modified type A 
coordination depicting the short barium-barium, separation, 
4.7 A, (dash line) and the small barium-oxygen-barium bond 
angle, 115°, (dot curve). 

The angle was calculated on the basis of a 2.8 A barium-oxygen separation and 
the 4.7 A barium-barium separation obtained from the radial distribution curve. 



oxygen separation, another is a 1.4 A boron-oxygen 
separation resulting from triangularly coordinated 
boron, and also present is a 1.5 A boron-oxygen sepa- 
ration arising from tetrahedrally coordinated boron 
[10]. The contributions of the two boron-oxygen 
separations to the area of this peak are indeterminate. 
However, compared with the contribution of the 1.6 
A silicon -oxygen separation, the contributions due to 
the boron-oxygen separations are minor because of 
the low atomic number of boron. With due consid- 
eration for these facts, an estimated area of 400 (elec- 
trons) 2 is attributed to the 1.6 A silicon-oxygen sepa- 
ration. From this area a coordination number of 4.3 
is calculated for the number of oxygen nearest neigh- 
bors about a silicon atom (see reference 3 for coordi- 
nation number calculations). This value is consist- 
ent with the predicted value, 4. 

The area of the 2.8 A peak has contributions from 
several interatomic separations as in the previous 
case. The primary contribution to this peak is from 
the barium-oxygen nearest neighbor separation. 
The area of the 2.8 A peak is estimated to be 2,500 
(electrons) 2 . If one were to assume 8-fold coordina- 
tion for the oxygen nearest neighbors about a barium 
atom, then the calculated area is approximately 
1,500 (electrons) 2 . 

The more important aspects of this research, the 
coordination change and estimate of percent of 
barium atoms in either coordination grouping, can be 
treated further. Continuing with the previous 
assumption that the barium atoms form a cubic 
array, an average barium separation of 5.58 A is 
calculated using the Warren and Pincus equation 
which states that the extent of immiscibility is 
inversely proportional to the cube of the modifier 
cation separation. The radial distribution function 
does not allow a determination of the relative 
amounts of barium atoms in type B and modified 
type A coordination but, as the contributing distances 
and average distance are known, the following 
calculation can be made to provide an estimate of the 
fraction of barium atoms in each coordination type. 

Letting a equal the fraction of barium atoms 
involved in type B coordination and (1 — a) the 
fraction involved in modified type A coordination, 
then 

a(6.8 A) 3 +(l-c0(4.7 A) 3 = (5.58 A) 3 
a-0.33. 

Since the total amount of barium oxide in the 
modifier-rich liquid is 25 mole percent, then approxi- 
mately 8.25 mole percent is associated with type B 
coordination and about 16.75 mole percent is involved 
in modified type A coordination. 

The ratio of mole percents of barium oxide (16.4%) 
to boric oxide (83.6%) in the binary barium borate 
system at the limit of miscibility of barium oxide 
in the modifier-rich liquid is approximately 0.20. 
For this system the barium atoms have the type B 
configuration [3]. It is interesting to note that the 
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ratio of mole percents of barium oxide (8.25%) to 
boric oxide (40%) for the corresponding modifier- 
rich liquid in the ternary barium borosilicate system 
is 0.21, when only the calculated amount of barium 
oxide involved in type B coordination is considered. 
The similarity of ratios suggests that the amount of 
barium atoms in type B coordination remains con- 
st nut as silica and barium oxide components are 
added, while the added barium atoms assume the 
modified type A coordination. 

Using the oxygen-volume method [1] the compo- 
sition of the maximum extent of miscibility of modi- 
fier oxide can be estimated. For barium borate in 
type B coordination, the relative mole percentages 
are: 16.1 BaO-83.9 B 2 3 [1]. For barium silicate in 
modified type A coordination, the relative mole 
percentages are: 32.2 BaO-67.8 Si0 2 . The combined 
compositions can be written as: 32.2 BaO (modified 
type A). 16.1 BaO (type B)-83.9 B 2 3 -67.8 Si0 2 or 
on a 100 percent basis as: 16.1 BaO (modified type 
A). 8.0 BaO (type B) • 42.0 B 2 3 -33.9 SiO,. This 
composition compares well with the experimental 
composition, 16.75 BaO (modified type A) -8.25 BaO 
(type B)-39 B 2 3 -36 Si0 2 , where the relative amounts 
of each coordination type were determined by the 
calculation of a, the fraction of barium atoms 
involved in type B coordination. 



4. Conclusions 

The radial distribution functions reveal barium- 
barium separations at approximately 4.7 and 6.8 A. 
The unexpected 4.7 A distance is due to barium 
atoms in a modified type A coordination with a 
calculated barium-oxygen-barium bond angle of 
about 115°. The 6.8 A distance is due to barium 
atoms in type B coordination and the barium- 
barium second nearest neighbor separation due to 
modified type A coordination. The new factor as 
compared to the average barium separation in the 
barium borate glass is the 4.7 A separation. As 
postulated by Levin and Block, the increase of the 
immiscibility range in the ternary barium boro- 
silicate system does occur by a coordination-change 
mechanism. However, the barium atoms are not 
involved solely in type A coordination. The factor 
which caused Levin and Block to assume 100 
percent change to type A coordination is the good 
agreement observed between the average barium- 
barium separation (5.58 A) determined by using the 
Warren and Pincus immiscibility equation and the 
geometrical separation (5.50 A) due to the type A 
configuration with a 180° barium-oxygen-barium 
bond angle. 

At the limit of miscibility in the binary barium 
borate system, the modifier-rich liquid contains 
16.4 mole percent barium oxide with the barium 
atoms involved primarily in type B coordination. 
The corresponding barium borosilicate modifier-rich 
liquid contains 25 mole percent barium oxide [1], 
Combining the results of the radial distribution 
study and the immiscibility equation indicates that 



approximately 16.75 mole percent barium oxide is 
involved in the modified type A configuration and 
8.25 mole percent barium oxide is associated with 
type B coordination. In the barium borate modifier- 
rich liquid as silica and barium oxide components 
are added, it may be considered that the amount of 
barium atoms in type B coordination remains con- 
stant, while the additional barium atoms assume a 
modified type A coordination. Using the modified 
type A coordination model and the oxygen-volume 
method of calculating the extent of miscibility of 
modifier oxide [1] results in a final limiting compo- 
sition of 16.1 BaO (modified type A) -8.0 BaO 
(type B)-42.0 B 2 3 -33.9 Si0 2 in mole percent. 

Each coordination type, modified A or B, is not 
necessarily associated with a particular network 
former. One should consider the structure with 
type B coordination as having modifier cations 
separated primarily by the network former as 
described earlier in figure 2 and related legend. Thus, 
the modifier cations are bonded to different oxygens. 
The structure with modified type A coordination 
should be considered as having modifier cations 
bonded primarily to the same oxygen atom. 

This study indicates the importance of determin- 
ing whether the phenomenon of coexistence of type 
B and modified type A coordinations is limited to the 
ternary system. Research on relevant binary 
systems is presently being undertaken. 



The authors are greatly indebted to A. Bienen stock 
for his valuable counsel in connection with many 
aspects of this study. In addition, thanks are due 
to H. M. Ondik of the National Bureau of Standards 
staff for developing a computer program for use in 
facilitating the calculations connected with this 
research. 
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